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Abstract. We characterize all linear operators on finite or infinite-dimensio- 
nal spaces of univariate real polynomials preserving the sets of elliptic, positive, 
and non-negative polynomials, respectively. This is done by means of Fischer- 
Fock dualities, Hankel forms, and convolutions with non-negative measures. 
We also establish higher-dimensional analogs of these results. In particular, 
our classification theorems solve the questions raised in |§] originating from 
entire function theory and the literature pertaining to Hilbert's 17th problem. 



1. Introduction 

The dynamics of zero sets of polynomials and transcendental entire functions 
under linear transformations is a central topic in geometric function theory with 
applications ranging from statistical mechanics and probability theory to combina- 
torics, analytic number theory, and matrix theory; see, e.g., [3]-[8] and the refer- 
ences therein. Describing linear operators that preserve non- vanishing properties is 
a basic question that often turns out to be quite difficult. For instance, the problem 
- going back to Laguerre and Polya-Schur |19| - of characterizing linear operators 
preserving univariate real polynomials with all real zeros remained unsolved until 
very recently [3]. Just as fundamental are the problems - originating from entire 
function theory [Til [121 [El [TTl [20] and the literature pertaining to Hilbert's 17th 
problem [ini [5T1 [55] - of classifying linear operators preserving the set of univariate 
real polynomials with no real zeros or the closely related sets of positive, respectively 
non-negative polynomials. Let us formulate these problems explicitly. 

Given d S No, where No is the additive semigroup of non-negative integers, let 
R^[x] = {peR[x]: deg(p) < d} and set 

POS = {pe R[x] : p{x) > 0, x e R}, POSd POS n Rd[x], 

SOS= i^J2pi{xf -.pi eR[xl 1 <i <r< ooj , SOSa = SOS DRaix], 

ELL = {p e M[x] : p{x) 7^ 0, a: e R}, ELL a = ELL n R<j[x]. 



It is well-known that SOS = POS := {p £ R[a;] : p{x) > 0, a; £ M}, ELL = 
POS^{-POS), and POS = \Ja>n{SOS + a). Polynomials in POS, SOS, ELL are 
called positive, non-negative or sos (sums of squares) , and elliptic, respectively. Let 
Zd be the set of all linear operators T : Rd[x] M\x\ and £ be the semigroup of all 
linear operators T : M.[x\ R[x]. For V e {POS, SOS, ELL} and deNo define 

£(y) - {T e £ : T{V) C V}, £d{V) = {T e £d : T{Vd) C V}. 

Problem 1. Characterize 2,d{,P0S), Zd{SOS), td{ELL) for each even d G No- 
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Problem 2. Characterize the semigroups £{POS), £{SOS), £{ELL). 

In physics it is useful to distinguish between between "soft" and "hard" theorems 
asserting the non-vanishing of partition functions in certain regionsQ By analogy 
with this terminology and the one used in the classification of linear preservers of 
real-rooted polynomials [31 H] , one may say that results pertaining to Problem [T] 
are "hard" or "algebraic" (bounded degree) while those for Problem [2] are "soft" or 
"transcendental" (unbounded degree). Special cases of Problems [THU were recently 
studied in f9], where one can also find a list of partial results known so far and 
corrections to some erroneous claims from the existing literature. 

In this paper we give complete solutions to Problems HHH In particular, this 
supersedes P| whose main results are now immediate corollaries of the classification 
theorems in iJ2J In [JS] we answer multivariate versions of Problems [THU while in ijH 
we establish some related results and discuss further directions. 

One of the tools that we use both in the univariate and the multivariate case 
is the following inner product structure on polynomial spaces. Recall (cf., e.g., [B] 
and the references therein) that the Fischer-Fock or Bargmann- Segal space Tn is 
the Hilbert space of holomorphic functions / on C" such that 



where z = (zi, . . . ,z„) G C", |z|^ = X]"=i kd^j ^^^d dV(^) is the volume element 
(Lebesgue measure) in C". The inner product in JF„ is given by 



Here a! — ai!---a„! for a = (ai,...,a„) G Nq while J2a^a2", respectively 
^^6q,z", is the Taylor expansion of /, respectively g, where z" = z"^ ••■z"". 
Clearly, the space of all polynomials in n variables with real coefficients M[zi, . . . , z„] 
is contained in and thus inherits a natural inner product structure. Given a 
cone K C R[zi, . . . , z„] we let K* C R[zi, . . . , z„] be its dual cone, i.e.. 



The Fischer-Fock space J-n was used by Dirac to define second quantization and 
its inner product has since been rediscovered in various contexts, for instance in 
number theory - where the corresponding norm is known as the Bombieri norm 
- and the theory of real homogeneous polynomials [IS] . In D-module theory and 
microlocal Fourier analysis one usually works with the inner product on J^n defined 
by (/(z), ,g(z)) = {f{iz),g{iz)). For these and further properties of ^„ (such as its 
Bergman-Szego-Aronszajn reproducing kernel) see, e.g., the references in [5]. 

Another main ingredient in our characterization theorems is the theory of posi- 
tive definite kernels, Hankel forms, and moment problems in one or several variables. 
There is a vast literature on these topics ranging from classical works on univariate 
and multivariate moment problems [U [TOl [TH [151 [23 [29] to recent contributions 
[21 [TBI [22j [23l [211 [26] and we refer to these for additional details. As usual, a k x k 
real symmetric matrix A is said to be positive semidefinite (respectively, positive 
definite) if x*Ax > for aU x G M'^ (respectively, x*Ax > for aU x G M*^ \ {0}). 



This dichotomy stems from "soft-core" vs. "hard-core" pair interactions in lattice-gas models 
and does not refer to the level of difficulty in proving such theorems but rather to the fact that 
in some sense "soft" theorems are constraint-free while "hard" theorems involve constraints of 
various kinds such as the maximum degree of a graph; see, e.g., [3] and the references therein. 





'0 



K* = {f &R[z,,...,z,,]:{f,g)>0,ge 



K}. 
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Let dJli denote the class of non-negative measures on R with all finite moments, 
that is, < oo, m G Nq. Hamburger's theorem [Uj asserts that a sequence 

of real numbers {am}m£No ^he moment sequence of a measure /i G , i.e., 

a™ = y t'^dfiit), m e No, 

if and only if the Hankel matrix (ai+j)™ ^g is positive semidefinite for any to g Nq. 
If such a ^ is unique the corresponding moment problem is said to be determined. 

Haviland's theorem [T5] is a multivariate analog of Hamburger's theorem that 
may be stated as follows (cf., e.g., [21 Chap. 6]). Let n G N and be the class of 
non-negative measures ^ on R" with all finite moments, that is, /|t"|d/x(t) < oo, 
a G NJJ, where t = (ii,.. .,i„) G M" and t" = t^^ ■■■t^" for a = (ai,. .. ,a„) G NJf. 
A sequence of real numbers {acjcgNjj satisfies 

flc, = y" t"d^(t), a G N;J, 

for some fi G Tl^ if and only if the following holds: the real-valued linear functional 
L on R[xi, . . . , Xn] determined by Llx") = a^, a G Ng , satisfies L{f) > for any 
non-negative polynomial / G M.[xi, . . . , Xn] (i.e., /(x) > 0, x G R"). 

Acknowledgement. I would like to thank the anonymous referee for interesting re- 
marks and useful suggestions. 

2. Univariate Polynomials 

2.1. Preliminaries. For convenience, we recall some elementary properties of the 
classes of linear operators under consideration. If T G £d and T(l) G ELL then the 
following are equivalent: (i) T G Zd{ELL); (ii) T G 2,d{P0S) or -T G £,d[POS); 
(iii) T G Z.d{SOS) or -T G £d(5'0S'). Similarly, if T G £ and r(l) G ELL then 
the following are equivalent: (i') T G £{ELL); (ii') T G 2{POS) or -T G i2(POS); 
(iii') T G £(S'OS') or -T G £(S'OS'). If one drops the assumption r(l) G ELL 
then (i) ^ (ii) => (iii) and (i') ^ (ii') => (iii'). Furthermore, if T G SllsOS) and 
T(l) = then T = 0. These are simple consequences of straightforward homotopy 
arguments (see, e.g.. Remark l4!2l below and [9]). 

2.2. Transcendental Characterizations. We start with the "soft" (unbounded 
degree) case and the three classification questions stated in Problem [2] Recall that 
any linear operator T d £ may be uniquely represented as a formal power series in 
D := with polynomial coefficients, i.e., 

oo 

T^Y.q,{x)D\ q,eR[x],ieNo, (2.1) 

1=0 

see, e.g., [S]. Given T as above we associate to it a one-parameter family of linear 
differential operators with constant coefficients {Tyj^gR defined as follows: 

oo 

Note that in general Ty has infinite order and Tj, = T, y G M, if T has constant 
coefficients. For each ?/ G K and to G No we also define a polynomial Pj,,m = Pt^ „ 
of degree at most to and an (to + 1) x (771 + 1) Hankcl matrix TCy.m ~ m t>y 

m 

Py,ra{x) = ^qt{y)x' G R[x], Uy^ra ^ ( (i + j) Iqj+j (?/ )) ™.^g . 
j=0 
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Remark 2.1. Clearly, py^rn{x) £ R[a;,?/] for each m g No. We call this polynomial 
in two variables the m-truncated symbol of the operator T S £. 

Our first theorem characterizes the semigroup £,{SOS) of linear operators pre- 
serving the cone SOS of non-negative/sos polynomials. 

Theorem 2.1. Let T ^ £, be as in (12. ip . The following assertions are equivalent: 

(1) T €Z{SOS); 

(2) Ty e £(505) for all y e R; 

(3) {py,27n, /) > for ally eR, m e No, / G S0S2rn, i-e., Py.m e 5052„; 

(4) tiy^m is positive semidefinite for all y d m G No; 

(5) For any s €R there exists fi^ G such that Jt^dnsit) G M[s], i G No, 
and i\qi{y) = Jt'^diJ,y{t), y G M, i G No- The family of operators {TyjygK is 
then given by the convolutions 

Ty{f){x) = j f{t + x)diiy{t), / G R[x], y G M. 

Proof. Clearly, (2) (1). If (1) holds then since SOS is invariant under shift 
operators e"-°, a G R, one gets e^^^Te"^ G £(505) hence 

5]g.(x-a)/«(x)>0 

1=0 

for any / G 505 and a, a; G R, which proves (2) (note that for each / G M.[x] the 
sum in the above inequality is finite). Next, if a, y G R and / G Rrf[a;], d G No, then 

r,(/)(a) = {PyM,e^''f). (2.2) 

i=Q 

Since / G SOSd ^ e"-°/ G SOSd, a G R, identity ([221) shows that if (2) holds then 
{Py.d, /) > for all y G R, 0? G No, / G SOSd- In particular, this is true for all even 
d, which proves (3). Conversely, if (3) holds then (2) follows from (|2.2p since any 
/ G SOS has even degree. Now, for m G No, y G R, and g{x) = J2™^o '^^^^ ^ RmN 
one has 

m 

{Pyar>i,9^) = X! (* + j)'*+j(2/)a»aj 

and thus (3) <^ (4) since any / G SOS is a sos. Finally, by Hamburger's theorem 
(see SjlJ the following holds for each y G R: y-y,m is positive semidefinite for all 
m G No if and only if there exists /ij, G 971]'" such that i\qi{y) = f t^d^y{t), i G Nq. 
Since all qi must be polynomials we conclude that (4) <^ (5), which completes the 
proof of the theorem. □ 

Remark 2.2. By Carleman's sufficient condition for the Hamburger moment prob- 
lem to be determined jTU], if y G R is such that X^i^i [(20''Z2i(y)]~^^^' = oo then the 
corresponding measure jjLy G 37lf in Theorem l2.1l (5) is unique. Recent progress on 
the determinateness of multivariate moment problems has been made in [221 [531 [51] . 

The following analog of Theorem 12.11 characterizes the semigroup Z{POS) of 
linear operators preserving the cone POS of positive polynomials. 

Theorem 2.2. Let T ^ 2, be as in (j2.ip . The following assertions are equivalent: 

(1) T G Z{POS); 

(2) Ty G Z{POS) for all y G R; 

(3) (Pa,2m, /) > for allyeR,m& No, / G POS2m; 

(4) tiy^rn is positive definite for all y G R, m G No; 

(5) det(7^2;,„0 > for ally eW, me No; 
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(6) For any s G M there exists /ig £ 9Jlf whose support supp(//s) is an infinite 
set such that Jf^d^sit) G ^[s], i & No, and i\qi[y) — Jt'-dfiy{t), y G R, 
i (z Nq. The family of operators {Tyj^gR is then given by the convolutions 

Tv{f){x) = jf(t + x)df,yit), f e R[x], y e M. 

Proof Recall that POS = Ua>o{SOS + a) and note that if T e Sl{POS) then 
go G POS. Using these observations and the fact that the cones POS, P0S2m, 
m € No, are invariant under shift operators e°-°, a € K, one then gets (1) <^ (2) 
<^ (3) (4) (6) as in the proof of Theorem O (see, e.g., [Ml Theorem 1.2] 
and Remark 2.3 in [2 Chap. 6] for the extra condition that the supports of the 
representing measures should be infinite sets, which is not required in Theorem 1 2. II 
(5)). The equivalence (4) <=> (5) is a classical result in matrix theory due to Jacobi 
and independently to Hurwitz, cf., e.g., [T51[2S]. □ 

By the remarks in H2.1\ Theorem 12.21 also yields a description of the semigroup 
£,{ELL) of linear operators preserving the set of elliptic polynomials. Problem [5] is 
therefore completely solved. 

2.3. Algebraic Characterizations. Since there are no positive, non-negative or 
elliptic polynomials of odd degree, for the "hard" (bounded degree) classification 
questions stated in Problem [1] it is enough to consider only linear operators T £ Zd 
with d = 2fc, fc G No. It is well-known that any such operator may be viewed as 
an element of the Weyl algebra ^i(M) of order at most d (cf., e.g., [Hj), that is, a 
linear ordinary differential operator with polynomial coefficients of the form 

d 

T = Y.q,{x)D\ q,eR[x],ieNo. (2.3) 

i=0 

By analogy with the "soft" case treated in H2.2I we associate to T a one-parameter 
family of linear differential operators of order at most d with constant coefficients 
{TyjygR, a one-parameter family of polynomials {Py,d}yeR in I^dN, and a one- 
parameter family of {d+1) x (d+ 1) Hankel matrices {'Hy.d}yeK defined as follows: 

d d 

Ty = ^(Z,(y)i?\ PyA^) = ^g,(y)a;^ G Rd[x\, Hy^ - {{i + j)\q,+j{y))l^^^. 

Remark 2.3. Given d G No it is clear that py d G M.[x, y\. This two variable polyno- 
mial is called the algebraic symbol of the operator T £ £,d- 

The finite degree versions of Theorem 12 . 1 1 and Theorem 12 . 21 given below describe 
the sets £,d{SOS) and £,d{POS), respectively. We omit their proofs since they are 
almost identical to those of Theorems I2.1H2.2I 

Theorem 2.3. For d = 2k and T £ £,d as in (|2.3p the following are equivalent: 

(1) TeSldiSOS); 

(2) Ty G £,diSOS) for all y G M; 

(3) {Py,d, /> > for all y G M, / G SOSd, i.e., Py^a e SOS*^; 

(4) Ti-yA is positive semidefinite for all y G M. 

Remark 2.4. In his work '27] on the truncated Hamburger moment problem M. Riesz 
showed that the Hankel matrix associated to a given finite sequence of real numbers 
is positive semidefinite if and only if the sequence satisfies certain equalities and 
inequalities involving moments of non-negative measures. Using this result one can 
add a fifth condition to Theorem 12.31 (equivalent to those already stated) which is 
similar to Theorem l2.1l (5). 
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Theorem 2.4. For d ^ 2k and T G £,d Cis in (|2.3p the following are equivalent: 

(1) T e WOS); 

(2) Ty € ZdiPOS) for all y G K; 

(3) {py,d, /) > for allyeR, fe POSa; 

(4) Hy^d is positive definite for all y G M; 

(5) det(Hy,„0 > /or an J/ G M, < m < d. 

Remark 2.5. In the case of finite order linear ordinary differential operators with 
constant coefficients Theorems I2.IH2.2I reduce to classical results of Hurwitz [17] 
and Remak [25], see also Polya-Szego's [20l Chap. VII]. 

Note that from Theorem 12.41 and the observations in §2.11 one also gets a char- 
acterization of the set Zd{ELL). Therefore, these results fully solve Problem[T] 

3. Multivariate Polynomials 

We will now establish multivariate analogs of the characterization theorems in 

..,a„) gN;J let 

.,Xn] ■■ deg^^(p) < tti, 1 < « < n}, 
. ,x„] : p[xi, . . . ,x„) > 0, (xi, . . . ,x„) G M"}, 
. ,x„] : p{xi, . . . ,x„) > 0, (xi, . . . ,a:„) G M"}, 



1 Given n G N and a = (ai 




^{pER[xi,. 


POS{n) 


= {peR[xi,. 


POS{n) 


= {peR[xi,. 


SOS{n) 


= < '^P^{Xl, 
U=l 


ELL{n) 


= {peR[xi,. 




^vnRa[xi, 



■ ,Xn)'^ ■ Pi &R[xi, . . . , Xn] , I < i < r < ooj , 

,x„] : p{xi, . . . ,x„) 7^ 0, (xi, . . . ,x„) G K"}, 

,x„], V £ {POS{n),'POS{n),SOS{n),ELL{n)}. 

These sets are the natural multivariate analogs of the ones defined in [JT] for n = 1 . 
Note that SOS{n) C POS{n) with strict inclusion in general. Let £,{n) be the 
space of all linear operators T : R[xi, . . . , Xn] — > IR[a;i, . . . , a;„]. As is well-known 
(cf., e.g., |B]), any T G may be uniquely represented as an (d priori infinite 

order) linear partial differential operator with polynomial coefficients: 

T=^<7„(x)a", g„ GE[a;i,...,x„], aGN^\ (3.1) 

where x = (xi, . . . , a;„), 9" = 9"^^ • • ■ for a = (ai, . . . , a„) G Nq and 9a;, = 
d/dxi, 1 < i < n. We then associate to T an n-parameter family of linear partial 
differential operators with constant coefficients {TyjygKri given by 

for any fixed y = {yi, . . . , yn) G M". Moreover, for all a = {ai, . . . , an) & Nq we 
define an n-parameter family of polynomials {Py.alygR" by 

/3<Q 

for each fixed y = [yi,. . . , yn) G M", where x'^ = x^^ ■ ■ ■ x^" for x = [xi,. . . , a;„), 
P — {Pi, . . . ,Pn) e Nq, and we use the standard (product) partial order on Ng: 
P < a provided that Pi < ai, 1 < i < n. Let us also introduce an n-parameter 
family of real-valued maps {7iy}ygRn on the additive semigroup Nq by setting 

Hy{a) = a\qa{y), a = (ai, . . . , a„) G Nq , 

for any fixed y = (yi, . . . , ?/„) G M", where as before a! = ail ■ ■ ■ a„!. 
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Remark 3.1. Clearly, py^a G K[a;i, . . . , Xn, yi, . ■ . , j/„] for each a G Ng . By analogy 
with the terminology used in the univariate case (cf. Remark 12. ip we call this 2n 
variable polynomial the a-truncated symbol of the operator T. 

Set 2a = {2ai, . . . , 2a„) £ Nq whenever a = (ai, . . . , q;„) G Nq . The following 
theorem characterizes all linear operators in £(n) that preserve the cone POS{n) 
of non-negative polynomials in n variables. 

Theorem 3.1. For T £ S,{n) as in (|3.ip the following assertions are equivalent: 

(1) T( POS{ n)) C POSj n); 

(2) TyiPOS{n)) C POS{n) for all y e R"; 

(3) (Py.2a,/) > for all y G K", a £ / £ P0^(n)2a, i.e., Py,2a 6 

(4) For any s = (si, . . . , s„) G M" t/iere exists fig G smc/i i/iat J t"(i/is(t) S 
]R[si,...,s„], Of e NJf, andHyia) = /t"(i/Lty (t), y e M", a e NJJ, w/iere 
t" = • • - i^" /ort = (ti, . . . G M". The family of operators {TyjygR^ 
is f/ien given by the convolutions 

ry(/)(x)- Jf{t + ^)dfiy{t), /eR[xi,...,a;„],yeR". 

Furthermore, each of these conditions implies: 

(5) r/ie map Wy : Nq — > M is a positive semidefinite kernel on Ng x Ng for each 
y G M", i.e., the Hankel matrix (Ti.y{a^ + a^)'j^_._-^ is positive semidefinite 
for any to G N and {a^ , . . . , a™} C Nq . 

Proof. For each a — (ai, . . . , q;„) G Nq the cone POS{n)a is invariant under shift 
operators e^'^, a — (ai, . . . ,a„) G K", where a • 9 = aid/dxi + • • • + Ondjdxn. 
Moreover, one has 

Ty(/)(a) - I0(y)d''f{^) = {Py.a,e''-^f), / G R„[a;i, . . . , x„]. 

/3<Q 

Note also that if / G POS{n) then Aeg^^{f) = mod 2, 1 < i < n. The equiva- 
lences (1) <J=> (2) <^ (3) follow readily from these facts and straightforward extensions 
of the arguments in the proof of Theorem 12. ll to the multivariate case. 

Next, since Ty(/)(0) = {py^p, f) if f ^ R/3[xi, . . . , x„], /3 G NJJ, Taylor's formula 
yields (4) (3). To prove the converse, for y G R" we let Ly : R[xi, . . . , Xn] R 
be the linear functional determined by Ly(x'^) = Hy{P), (3 G Ng. It follows that 

Ly{f) = (Py./3, /), / e ^p[xi, . . . , a;„], /? G W^. 

In particular, if (3) holds then Ly{f) > for aU / G 'POS{n)2a, a G NJf. From 
Haviland's theorem (cf. ^ and the fact that 7iy(/3) is a polynomial in y for each 
P G we deduce that (3) ^ (4). 

Let now to G N, {«!,... , a'"} C c, G R, 1 < i < TO, and set 

m 

5(x) = ^ Qx"' G R[a;i , . . . , a;„] . 
1=1 

If a G Ng is such that a* -I- < a, 1 < z, j < to, one has the identity 

m 

(Py,„,<?2) = ^ Hy(a' +a^)c,Cj. (3.2) 

Since G POS{n) this proves that (3) =^ (5). □ 

The next theorem gives a characterization of all linear operators in £{n) that 
preserve the cone POS{n) of positive polynomials in n variables. 
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Theorem 3.2. For T G 2,{n) as in p.ip the following assertions are equivalent: 

(1) T{POS{n)) C POS{n); 

(2) Ty{POS{n)) C POS'(n) for all y e M"; 

(3) (Py.2a, /) > /or aZ/ y e R", a G N^, / e P05(n)2„; 

(4) For any s — (si, . . . , s„) G M" t/iere exists G whose support is not 
contained in a proper real algebraic variet'^ such that for all a G Ng one has 
J t"d^s(t) G R[si, . . . , s„] and Hy(a) = / t"d^y(t), y G M". T/ie /amiZ?/ of 
operators {TyjygK" *s t/ien given by the convolutions 

ry(/)(x)= y/(t + x)d/.y(t), /GM[xi,...,x„],yGlR". 

Furthermore, each of these conditions implies: 

(5) r/ie map Hy : Nq ^ E is a positive definite kernel on Ng x Nq /or eac/i 
y G M", i.e., f/ie Hankel matrix {Hyid^ + a-'))™^]^ is positive definite for 
any to G N and {a^, . . . , a™} C Ng . 

Proof. Since POS{n)2a is invariant under shift operators e'*'^, a G R", the ar- 
guments in the proof of Theorem 13.11 carry over mutatis mutandis to the present 
situation so we wiU not repeat them. The only novelty concerns the extra condition 
in (4) above on the support of the measures fig (this condition does not appear in 
Theorem 13. II f4)) for which we refer to e.g. Exercise 3.14 in [51 Chap. 6]. □ 

Remark 3.2. Elementary homotopy arguments show that ELL{n) = POS(n) U 
{-POS{n)) and if T G £(n) then T{ELL{n)) C ELL{n) if and only if either 
T[POS[n)) C POS{n) or T{POS{n)) C -POS{n). Therefore, Theorem [O also 
yields a description of all linear operators in £(n) preserving the set ELL{n) of 
elliptic polynomials in n variables. 

The above theorems solve the multivariate analogs of the "soft" (transcenden- 
tal/unbounded degree) classification questions stated in Problem [21 It is clear that 
one gets "hard" (algebraic/bounded degree) theorems in similar fashion. More pre- 
cisely, given a G Ng and a linear operator T : R2Q[a;i, . . . , a;„] — > R[a;i, . . . , Xn\ it 
follows from the proofs of Theorems I3.1H3.2I that 

T{P0S{n)2o.) C TOSin) ^ Ty(POS{n)2a) C PaS(n) for aU y G R" 

^ (Py,2a, /) > for all y G R", / G P0^(n)2a 

(i.e., Py'2a G POS{n)*^ for aU y G R") and 

T{P0S{n)2a) C POS{n) ^ Ty{POS{n)2c,) C POS{n) for aU y G R" 

^ (Py,2a, /) > for aU y G R", / G P0S{n)2a- 

These equivalences answer the multivariate analogs of the classification questions 
raised in Problem [Tl 

4. Related Results and Further Directions 

1. It is interesting to compare the classification theorems established here with 
those obtained in (31 [H [5] for linear operators preserving real-rooted polynomials 
or their multivariate analogs. A common feature is the use of appropriately defined 
operator symbols, see Remarks 12.11 12.31 and 13. II However, the truncated symbols 
introduced here are employed in a different way than the symbols used in [3l [4] . 
Note for instance that if T G £,{POS) is as in (|2.ip then for each to G Ng its 2to- 
truncated symbol Py^2m{x) is a positive polynomial in variables x, y (since all even 
degree truncations of the exponential series are positive) but the converse statement 



^In other words, supp(/Xs) is not a subset of / -"^(O) for some / S K[xi, . . . ,Xn\ \ {0}. 
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is not true, as one can see from e.g. Corollary [ITT] (b) below. By contrast, the main 
results in ^ essentially assert that a linear operator preserves the set of univariate 
real-rooted polynomials if and only if its symbol is a "real-rooted" polynomial in 
two variables (in some appropriate sense). 

We should also note that there are plenty of linear differential operators of pos- 
itive order in the Weyl algebra ^i(M) which preserve real-rootedness (these are 
characterized in [l^iBj) whereas no such operators preserve positive, non- negative 
or elliptic polynomials (cf. Corollarv l4.ip . 

The equivalent conditions in our main results successively describe linear op- 
erators preserving non-negativity by means of linear differential operators with 
constant coefficients, Fischer-Fock dual cones, positive semidefinite Hankel forms, 
and convolution operators induced by non-negative measures (Theorems 12.11 and 
I3.ip . The Fischer-Fock product was also used in Theorem 1.11 of [B] which asserts 
that an operator T G Ai (M) preserves real-rootedness if and only if its Fischer-Fock 
adjoint T* has the same property. Note that this duality fails for linear operators 
preserving positive, non-negative or elliptic polynomials since multiplication by a 
positive polynomial clearly preserves the sets SOS, POS, ELL but differential 
operators of positive order in ^i(R) do not. Nevertheless, Theorem 12.11 (3) and 
Theorem 13.11 (3) show that appropriate Fischer-Fock dualities (in one or several 
variables) do hold in the present context as well. 

2. The classification theorems in i j2.2H2.3l vield in particular the following improve- 
ments of the main results in [9J (Theorem A, Corollary, and Theorem B in op. cit.). 

Corollary 4.1. In the above notations one has: 

(a) IfT E £ is an element of the Weyl algebra Ai{M) of order d > 1 then for any 
even integer i > d one has T(yi)\V ^ whenever V £ {POS, SOS, ELL}; 

(b) There are no linear ordinary differential operators with polynomial coeffi- 
cients of finite positive order in £{POS) U 2{SOS) U £,{ELL); 

(c) // T G £ \ {0} is a linear ordinary differential operator with constant co- 
efficients then T e £{POS) <^ T e £.{SOS) ^ T e 2{ELL) ^ T is a 
convolution operator of the form 

Tif){x) = Jfit + x)d^i{t), feRix], 
for some ii € 971^ . 

Proof. To prove (a) it is enough to show that T ^ £^(F) for the smallest even 
integer i > d and V — SOS (the arguments for V — POS and V = ELL are 
similar). Let 

d 

T = Y, q^{x)D\ q, e R[x], I e No, qd ^ 0, 

and y G R be such that qd(y) 7^ 0. If d is odd set q^j^xiy) = 0. Then the determinant 
of the lower right 2x2 principal submatrix of Ti-yA+i is —qdiv)^ so TLy^d+i is not 
positive semidefinite and by Theorem l2.3i r ^ £,d+i{SOS). If d is even set qd+i{y) = 
qd+2{y) = and suppose that T G Zd{SOS) (if T (f. Zd{SOS) there is nothing to 
prove since SOSd C SOSd+2)- By Theorem 12. 31 Ti.^ ^ is positive semidefinite hence 
<ld{y) > 0. The determinant of the lower right 3x3 principal submatrix of 'Hy_d+2 
equals —qdiu)^ and thus TCy^d+2 cannot be positive semidefinite. By Theorem 12.31 
again we conclude that T ^ £.d+2{SOS). 

Clearly, (a) (b) while (c) follows from Theorem 12.11 and the fact that if T G 
£{SOS) and r(l) = then T = (cf. gSI]). □ 
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Remark 4.1. Corollary 14. II (a) asserts that a linear ordinary differential operator of 
order d > 1 with polynomial coefficients cannot preserve positivity, non-negativity, 
or ellipticity when acting on Mf[x], where £ is any even integer greater than d. In 
particular, this implies part (b) stating that there arc no linear ordinary differential 
operators of finite positive order (acting on R[x]) that preserve positivity, non- 
negativity, or ellipticity. Finally, part (c) asserts that a non-trivial linear ordinary 
differential operator with constant coefficients (acting on M.[x]) preserves each of 
the sets POS, SOS, ELL provided that it preserves at least one of them and any 
such operator is a convolution with a non- negative measure with all finite moments. 



3. From Theorems I3.1H3.2I we deduce that any linear partial differential operator 
with constant coefficients that preserves non-negativity or positivity is actually a 
convolution operator induced by a non-negative measure with all finite moments. 
More precisely, we have the following multivariate analog of Corollary 14.11 (c) . 

Corollary 4.2. Let T e ^{n) he a linear partial differential operator with constant 
coefficients. 

(i) The following assertions are equivalent: 

(a) T(POS{n)) C POSin); 

(b) There exists /i G SOT^ such that T may be represented as the convolution 

T(/)(x) = y"/(t + x)dA*(t), f eR[xi,...,Xn]. 

(ii) The following assertions are also equivalent: 

(c) T{POS[n)) C POS{n); 

(d) There exists /i G 971^ whose support is not contained in a proper real 
algebraic variety such that T may be represented as the convolution 

T(/)(x) = y"/(t-Kx)d/.(t), f eR[x,,...,Xn]. 

A similar corollary follows for diagonal operators in the standard monomial basis 
of K[a;i, . . . ,Xn], i.e., operators T £ £,{n) given by T{x^) = A/jx^, where A/j G K, 
/? G Nq. It is easy to see that such an operator may be uniquely represented as 

r = ^ a„x"5°, a„ G M, a G NJ,\ 

and that Xp = J^aK/sWa'^a for (3 G Ng, where as usual (/?)„ = (/3i)ai ■ • ■ {Pn)a„ 
and {m)k = m{m — 1) ■ • • (m — A: -I- 1) whenever fc, to G Nq with k < m. Elementary 
computations then yield the following implication: 

a!aa = y"t"d^(t), a G NJ,' =^ A;, = j [t + if dpi{t) ^ ^ s'^di/(s), /3 G NJJ, 

where 1 denotes the all ones vector and s = t -I- 1 (note in particular that for /i, v 
as above one has /i G 9Jt^ ^ G 9Jl^). Hence, if the condition in the left-hand side 
of the above implication holds we deduce that the diagonal operator T is given by 

T(/)(x) = j f[sx)dy{s), f G M[a;i, . . . , x„], 

where sx = (sia;i, . . . , s„a;„). The next statement is therefore an immediate conse- 
quence of these arguments and Theorems I3.1H3.2I 

Corollary 4.3. Let T G S,{n) be a diagonal operator given by T{x^) = X/sx^ , 
A/3 G R, /5 G NJf. 

(i) The following assertions are equivalent: 
(a) T(POS{n)) C POSin); 
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(b) There exists v G such that T may he represented as 

T(/)(x) = J /(sx)di.(s), / e M[xi, . . . , Xn]. 

(ii) The following assertions are also equivalent: 

(c) T{POS{n)) C POS{n); 

(d) There exists v G whose support is not contained in a proper real 
algebraic variety such that T may be represented as 

T(/)(x) - J /(sx)di.(s), / e M[xi, . . . , Xn]. 

Remark 4.2. Let T e be either a partial differential operator with constant 

coefficients or a diagonal operator. Note that r(l) is then a constant polynomial 
and suppose that T(POS{n)) C 'POS{n). If r(l) = it follows from CoroUaryHJ 
(i) (respectively, Corollary 14.31 (i)) that the the non-negative measure pi (respec- 
tively, v) has total mass zero and then T = 0. Therefore, if T ^ then T(l) > 
and so if / e POS{n) one gets T{f) = T{f - a) + aT{l) e POS{n), where 
a > is such that f — a e POS{n) (clearly, such a exists for any / G POS{n)). 
We conclude that if T € £(n) \ {0} is as in Corollary [321 or Corollary then 
T(POS{n)) C POSin) ^ T{POS{n)) C POS{n). 

4. If n > 2 one has in general SOS{n) C POS{n) and it is natural to ask for a 
characterization of the set {T e £,{n) : T{SOS{n)) C POS{n)} (obviously, this set 
contains the semigroup of all T G ^{n) such that T{POS{n)) C POS{n) that we 
already described in Theorem 13. ip . The following theorem provides an answer to 
this question. 

Theorem 4.4. For n > 2 and T G S,{n) as in the following are equivalent: 

(1) T{SOS{n)) C TOSj n); 

(2) Ty{SOS{n)) C POS{n) for all y G K"; 

(3) (Py,2a,/) > for all y G M", a G NJf, / G SOS{n)2o., i.e., py,2a 6 
SOS{n)l^; 

(4) The map TCy : Nq M. is a positive semidefinite kernel on Ng x Ng /or each 
y G M", i.e., the Hankel matrix (7iy(a* -t- a^))™._j^ is positive semidefinite 
for any to G N and {a-'^, . . . , a™} C Nq . 

Proof. Since SOS{n) = Uae'N'i; SOS{n)2a and each cone SOS{n)2a is invariant 
under shift operators e"*'^, a G K", the imphcations (1) (2) ^ (3) ^ (4) follow 
as in the proof of Thcorcm l3.1l To show that (4) =^ (3) note that / G SOS{n)2a if 
and only if / = J2]=i 5j some r G N and gj G Ma[a;i, ■ • ■ , x„], I < j < r. Given 
a G Nq let m = TO(a) be the cardinality of {(3 G Nq '■ P < a} (which is clearly 
finite) and order this set lexicographically as {a^, . . . , a™}. Any g G KQ,[a;i, . . . , Xn] 
may then be written as 

■m 

f;(x) = ^ Cix" , Ci G K, 1 < i < TO, 

i=l 

so the desired conclusion follows from (|3.2p . □ 

In view of Theorem 14.41 and Corollary 14.21 one may ask whether the following 
analog of Hilbert's 17th problem holds, cf. [7', Problem 17]. 

Question. Is it true that for any p G POS{n) there exist f G SOS{n) and a 
linear partial differential operator with constant coefficients T G 2,{n) such that 
T{SOS{n)) C 'POS{n) and T{f) = p? 
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